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ABSTRACT: Huoyan goose is a Chinese local breed famous for its higher laying performance, but the problems of variety 
degeneration have emerged recently, especially a decrease in the number of eggs laid. In order to better understand the molecular 
mechanism that underlies egg laying in Huoyan geese, gene profiles in the pituitary gland of Huoyan geese taken during the laying 
period and ceased period were investigated using the suppression subtractive hybridization (SSH) method. Total RNA was extracted 
from pituitary glands of ceased period and laying period geese. The cDNA in the pituitary glands of ceased geese was subtracted from 
the cDNA in the pituitary glands of laying geese (forward subtraction); the reverse subtraction was also performed. After sequencing and 
annotation, a total of 30 and 24 up and down-regulated genes were obtained from the forward and reverse SSH libraries, respectively. 
These genes mostly related to biosynthetic process, cellular nitrogen compound metabolic process, transport, cell differentiation, cellular 
protein modification process, signal transduction, small molecule metabolic process. Furthermore, eleven genes were selected for further 
analyses by quantitative real-time PCR (qRT-PCR). The qRT-PCR results for the most part were consistent with the SSH results. Among 
these genes, Synaptotagmin-1 (SYTl) and Stathmin-2 (STMN2) were substantially over-expressed in laying period compared to ceased 
period. These results could serve as an important reference for elucidating the molecular mechanism of higher laying performance in 
Huoyan geese. (Key Words: Huoyan Goose, Pituitary Gland, Laying, Ceased, Suppression Subtractive Hybridization, Quantitative 
Real-time PCR) 



INTRODUCTION 

The Huoyan goose was listed as one of the seventy- 
eights nationally protected domestic animals by the Chinese 
goverrmient in the year of 2000. It was used for both meat 
and egg production, and was famous for its higher laying 
performance. But breed selection and conservation work for 
this famous Chinese local breed is relatively weak; the 
problems of variety degeneration have been very prominent, 
especially the number of laying eggs decreased. 

In poultry, the reproductive endocrine system and the 
reproductive activity are strictly controlled by the 
hypothalamic-pituitary-gonadal axis (Padmanabhan et al., 
2002). The laying performance is determined by the number 
of folhcles destined for ovulation and the capacity of the 
oviduct to transform the ova into a hardshelled egg. The 
egg-production process, from follicle maturation to 
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ovulation, is strictly controlled by hormones and other 
physiological conditions, including energy metabolism and 
the regulating of apoptosis of oocyte and follicle-associated 
cells. The activation and maintenance of normal follicular 
functions are dependent on gonadotropins including 
luterinizing hormone (LH) and prolactin (PRL) secreted by 
the pituitary, which in turn is regulated by the hypothalamic 
gonadotropin releasing hormones (GnRH) (Kuo et al., 
2005). Minute differences in hypothalamic or pituitary 
functions, might affect reproductive activities, including 
folhculogenesis, ovulation, oviposition, and incubation 
behavior (Etches et al., 1984). 

The goose is a seasonal breeder (laying during the short 
photoperiodic season) with a low egg production rate. The 
pituitary LH and PRL together with other hormones, 
function to synchronize or regulate the reproductive 
seasonality in both long and short day breeding birds (Shi et 
al., 2007). Consequently, the pituitary gland is ideal tissues 
to elucidate the molecular mechanism associated with egg 
production, since genes may be transcribed differently 
under different physiological conditions. Information of 
gene profiles in pituitary gland at different reproduction 
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period is essential for studying the molecular genetic 
mechanisms of egg laying and their application to 
improving the laying performance. 

Suppression subtractive hybridization (SSH) is a 
method to get a gene expression profile specifically 
designed for comparing gene expression at various 
developmental stages or in altered physiological conditions 
(Diatchenko et al., 1996). This powerful technique enables 
comparing mRNA of two populations and then obtaining 
clones of genes that are differentially expressed in one 
population but not in another. First, both mRNA 
populations are converted into cDNA, which the cDNA that 
contains differentially expressed transcripts as tester, and 
the reference cDNA as driver. Tester and driver cDNAs are 
hybridized and then the hybrid common sequences are 
removed. Consequently, the remaining unhybridized 
cDNAs represent genes that are differentially expressed in 
the tester yet absent from the driver mRNA. This allows 
isolation of genes that are specifically expressed under 
diverse sets of conditions. Based on this technology, nine 
transcripts related to high egg production in the chicken 
hypothalamus/pituitary gland were identified (Shiue et al., 
2006). As for the goose, 18 known and 8 unknown 
differentially expressed genes in the ovaries of egg-laying 
and pre-laying Zi geese were identified by this method 
(Kang et al., 2009). In addition, it was also used to detect 
differential expression of genes in the livers of laying and 
pre-laying geese (Ding et al., 2007). Although the 
differential expression profile of genes in the pituitary gland 
has been studied in laying geese compared with pre-laying 
geese (Yen et al., 2006), but to the best of our knowledge, 
information of Huoyan geese in laying period and ceased 
period has not been investigated so far. 

In this study, gene expression profiles in the pituitary 
gland of laying period and ceased period Huoyan geese 
were investigated using SSH coupled with qRT-PCR 
methods. Combined results will be useful in elucidating the 
molecular mechanisms that underlie egg laying and 
improving the laying performance of Huoyan geese. 

MATERIALS AND METHODS 

Geese and tissue collection 

The animal protocol used in this experiment was 
approved by the Animal Welfare Committee of Liaoning 
Province, China. The Huoyan geese were selected from 
Panjin Jiyuan goose breeding farm and raised according to 
the program used at this farm. During the experiment, geese 
were fed ad libitum with rice grain and were supplemented 
with green grass or water plants whenever possible. Feed 
was given during the daytime when the geese were released 
into an open area outside the house. Twenty female geese 



were killed by exsanguinations in January to obtain 
pituitary samples of ceased period geese. Another twenty 
female geese were killed in June to obtain pituitary samples 
of laying period geese. All pituitary samples were quickly 
dissected, frozen in liquid nitrogen, and stored at -80°C 
until total RNA extraction. 

Total RNA isolation and reverse transcription 

Total RNA was prepared by TRIzol reagent (Invitrogen 
Corporation, Carlsbad, CA, USA) according to the 
instruction of the manufacturer. The total RNA samples 
from ceased period and laying period pituitary (n = 20, for 
each) were pooled separately. The first- strand cDNA and 
ds-cDNA were synthesized using SMARTer™ PCR cDNA 
Synthesis Kit (Clontech, Mountain View, CA, USA) and 
later, the ds-cDNA was purified with QIAquick PCR 
Purification Kit (Qiagen, Hilden, Germany). 

Construction of the suppression subtractive hybridization 
cDNA library 

The cDNA libraries were constructed by SSH using a 
PCR-Select™ cDNA Subtraction Kit (Clontech) following 
manufacturer protocol. A forward SSH library was 
constructed to isolate the up-regulated genes of the laying 
period. It was used to identify clones in which the laying 
period cDNA was used as the tester and the ceased period 
cDNA as the driver. In contrast, a reverse SSH library was 
constructed to detect the down-regulated genes of the laying 
period. It was used to identify clones in which the ceased 
period cDNA was used as the tester and the laying period 
cDNA as the driver. After digestion with Rsal, the tester 
cDNA was divided into two subpopulations, and then each 
of which was ligated to a unique adaptor at 16°C for 8 h. 
The adaptor 1 -ligated or adaptor 2R-ligated tester cDNA 
was separately hybridized at 68°C for 8 h with an excess of 
driver cDNA after denaturation at 98°C for 90 s. The two 
hybridized samples were then mixed together without 
denaturation and hybridized once more at 68 °C for 12 h 
with excess of denatured driver cDNA. The resulting 
mixture was added with 200 yl dilution buffer and 
amplified by two rounds of suppression PCR. The primary 
PCR was performed using PCR primers in the kit. After the 
reaction mixtures were incubated at 75 °C for 5 min to 
extend the adaptors, PCR was performed at 94°C for 10 s, 
66°C for 30 s, and 72°C for 90 s for 27 cycles in a reaction 
volume of 25 |il. The PCR product was then diluted to 10- 
fold, from which 10 diluted liquid was used as template in 
the subsequently nested PCR. This time, the reaction was 
performed for 12 cycles of 94°C for 10 s, 68°C for 30 s, and 
72°C for 90 s, using nested PCR primer 1 and primer 2R. 
The PCR products obtained from the forward and reverse 
subtractions were separately cloned into the pMD-18-T 
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vector (TaKaRa, Dalian, China), and transformed into 
Escherichia coU DH-5a to estabUsh different cDNA 
libraries. 

DNA sequencing and bioinformatics analysis 

Positive clones were selected by conventional blue- 
white screening. White clones were randomly selected from 
both forward and reverse libraries. The positive clones were 
fiirther validated by colony PCR using nested PCR primers 
provided in the kit. The resulting products were subjected to 
the direct sequencing with Ml 3 primers. The vector 
sequences were removed through a Perl script and checked 
through VecScreen (http://www.ncbi.nlm.nih.gov/ 
VecScreen/VecScreen.html). Sequence analysis and 
homology comparisons were carried out using the BLASTn 
program of NCBI (http://blast.ncbi.nlm.nih.gov/Blast.cgi). 
Gene Ontology (GO) annotations were performed online 
from GORetriever (http://www.agbase.msstate.edu/). If no 
annotation was returned, GOanna was used to retrieve GO 
annotations assigned on the basis of sequence similarities. 
The resulting annotations were summarized according to 
the generic GOSUm set using GOSlimViewer (McCarthy et 
al., 2007). 

Quantitative RT-PCR 

To validate the differentially expressed genes detected 
by the SSH approach, eleven ESTs (expressed sequence 
tags) which six from forward library and five from reverse 
library were selected for qRT-PCR analysis. The qRT-PCR 
primers were designed using Primer 3.0 
(http://frodo.wi.mit.edu/primer3). All the information of the 
primers is listed in Table 1. Total RNA was extracted using 
TRIzol Reagent (Invitrogen Corporation, Carlsbad, CA, 
USA) according to the manufacturer's instructions. The 
concentration and purity of the RNA were measured using a 
spectrophotometer. Two micrograms of total RNA was 
reverse transcribed using PrimerScript RT reagent Kit 
(TaKaRa, Dalian, China). Real-time PCR was carried out 
on LightCycler® 480 II real-time PCR system (Roche). 



Each 25 )il reaction volume contained 1 )il 10 |iM (each) 
forward and reverse primers, 12.5 (xl 2xSYBR® Premix Ex 
Taq™ II (Takara, Dahan, China), and 2 |il cDNA products, 
and the final volume was adjusted using PCR-water The 
following PCR program was used for amplification: 15 min 
at 95°C, 40 cycles of denaturation at 95°C for 10 s and 
annealing and extension at 60°C for 30 s. Relative gene 
expression was calculated using the 2"'^'^'^' method, using the 
ceased period group as calibrator (relative expression =1) 
and P-actin as an internal reference gene. 

Statistical analysis was carried out by student's t tests 
with SPSS 16.0 for Windows (SPSS Inc.) to evaluate 
whether the means were significantly different (p<0.05). 

RESULTS 

Differential screening and EST sequencing 

Based on the results of the differential screening, 112 
cDNA clones from the forward library and 110 cDNA 
clones from the reverse library were randomly picked and 
sequenced. After trimming, 103 and 104 high quaUty ESTs 
from the forward and reverse library, respectively, were 
obtained, all the processed ESTs were submitted to NCBI 
for BLASTn. 

Analysis of BLAST results revealed that in forward 
library, as shown in Table 2, 30 sequences have high degree 
of similarity with known or predicted genes in GenBank, 
such as P21 protein (Cdc42/Rac) -activated kinase 2, 
transcript variant 2 (PAK2), Stathmin-2 (STMN2), 
Synaptotagmin-l(SYTl), Sprouty-related, EVHl domain 
containing 1 (SPREDl), Protocadherin 9 (PCDH9) and 
Vesicle-associated membrane protein-associated protein A 
(VAPA). In reverse library, as shown in Table 3, 24 
sequences have high degree of similarity with genes like 
Vesicle-associated membrane protein 4 (VAMP4), 
Calmodulin binding transcription activator 1 (CAMTAl), 
SPARC related modular calcium binding 2 (SMOC2), 
Cyclin-dependent kinase 5, regulatory subunit 1 (CDK5R1) 



Table 1. Primers used for the quantitative real-time PCR analysis 


Putative gene symbols 


Forward primer sequence (5 '-3') 


Reverse primer sequence (5'-3') 


Amplicon size (bp) 


CAMTAl 


TTTCATATCTCCTGCTGGAAG 


AGCAATATACCATGACACCC 


111 


ENPP2 


CACTGGGTGGATGATGAC 


TGCACGGAAACCATCAAC 


102 


SM0C2 


AGTGTGTAAATTGCTGTTGAC 


GTTGCATTCTCATTGCATGAT 


128 


VAMP4 


ACCATAAATAATGTCTTCGGGA 


CAGAAGAGAAACGTGCCACTA 


134 


CDK5R1 


ATGACACAGTCTGTCAGG 


GTTGACCAAATTGGTGGG 


115 


PAK2 


TCACACAGCAAGCACTGA 


TTAGTGTACTTAGGAGACGGT 


105 


PCDH9 


GAGAGTCAAGTAGAGGTCCATA 


TCTGATTTCTGCATTGCTTCC 


105 


SPREDl 


ACAAGAAAGACGTTACCACT 


GCTGCTACTCTGGTTGATT 


109 


STMN2 


CCGAGGTAAGAGGCAACAAATA 


TGGCTGTTTGGCATCATTTAG 


128 


SYTl 


TTCAGTCTGGCAGTGTAGG 


TCAGTTGTTCATTCAGGATAGC 


101 


VAPA 


TGCAACTGACACTTGCTAT 


GCTTTGCCTTTAATGATCTCG 


111 


ACTB 


CCTGGCACCTAGCACAATGA 


TGGGTGTTGGTAACAGTCCG 


194 



924 
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Table 2. Identities of differentially expressed eDNA clones from the SSH forward library 



iN^^oi oi^y\o 1 izenc name 


Accession nimiber' 


Species 


li-\ alue 


Contactin associated protein-like 2 (CNTNAP2) 


NM_001193337.1 


Taeniopygia guttata 


3.00E-23 


P2 1 protein (Cdc42/Rac)-activated Isinase 2, transcript variant 2 (PAK2) 


XM_422671.2 


Gallus gallus 


1 .OOE-92 


Phytanoyl-CoA 2-hydroxylase interacting protein-like (PHYHIPL) 


NM_001199504.1 


Gallus gallus 


l.OOE-97 


Solute carrier family 5 (sodium/glucose cotransporter), member 11 (SLC5A11), 


XM_414862.3 


Gallus gallus 


0.000006 


Serine/arginine-rich splicing factor 3 (SRSF3), 


JN M_OU 1 1 y 5 5 5 4 . 1 


Gallus gallus 


2.UUb-171 


Legumain (LGMN) 


AM_42i J25.J 


Gallus gallus 


1 AAHj qa 


otatnmm-z-iiKe (oiiviiNZj 


AM_UU JZUj 1 jZ. i 


Meleagris gallopavo 


O AATJ 


0-linked N-acetylglucosamine (GlcNAc) transferase 


jNM_0U10UoJ 17. z 


Meleagris gallopavo 


0 


(UDP-N-acetylglucosamineipolypeptide-N-acetylglucosaminyl transferase) (OGT) 








Glutamate-ammonia ligase (GLUL) 


■\T\/f On^AQI 1 

iNivi_zui>'fyj. 1 


Gallus gallus 


o.uuri-iuy 


Cytochrome c oxidase subunit I (COXl) gene, mitochondrial 


kju i /yuuz. 1 


Anser 


O.UUri-l /U 


Cytochrome c oxidase subunit III (COX3) gene, mitochondrial 


riTT17QAlA 1 


Anser 


A 


oynaptotagmm- i -luce a i i i ^ 


VM nn'?'7n9in9 i 

AiVl_UU JZUZ 1 UZ. i 


Melcagns gallopavo 




Ribosomal protein, large, PO (RPLPO) 


NM_2049o7.1 


Gallus gallus 


l.OOE-171 


Poly(A) binding protein interacting protein 2 (PAIP2) 


NM_00 1 007 8 3 2 . 1 


Gallus gallus 


8.00E-I09 


Minichromosome maintenance complex component 6 (MCM6) 


MA /t AA 1 AA(^CT7 1 

JN JVl_tXJ 1 (XJ032 / . 1 


Gallus gallus 


n AATir 1 CA 


Sprouty-related, EVHl domain containing 1 (SPREDl), 


NM_152594.2 


Homo sapiens 


2.00E-49 


Protocadherin 20 (PCDH20) 


XM_003980439.1 




8.00E-14 


Protocadherin 9 (PCDH9) 


XM_()()3640560. 1 


Gallus gallus 


0 


Vesicle-associated membrane protein-associated protein A-like (VAPA) 


XMJ)03204959.1 


Meleagris gallopavo 


0 


Solute carrier family 4, sodium bicarbonate cotransporter, member 4, (SLC4A4) 


XMJ)03641192.1 


Meleagris gallopavo 


l.OOE-161 


ATP synthase, H-H transporting, mitochondrial Fl complex, O subunit (ATP50) 


XM_416717.3 


Meleagris gallopavo 


0 


Four and a half LIM domains 2 (FHL2), 


XM_4I6924.3 


Meleagris gallopavo 


l.OOE-112 


Major facilitator superfamily domain containing 2A (MFSD2A) 


XM_4I7826.3 


Meleagris gallopavo 


0.0001 


SDAl domain containing 1 (SDADl) 


XM_420597.3 


Meleagris gallopavo 


0.0005 


Anser albifirons mitochondrion, complete genome 


AF363031.1 


Anser 


7.00E-170 


anser clone goosePiSSHFmixIBOS pituitary gland-expressed 


DQ836037.1 


Anser 


8.00E-55 


unknown gene 8 mRNA sequence 








Voucher BISE-Aves24 cytochrome oxidase subunit 1 (COT) gene.mitochondrial 


GU57 1728.1 


Anser anser 


6.00E-165 


Meleagris gallopavo myelin proteolipid protein-like 


X]yi_003205647.1 


Meleagris gallopavo 


0 


Chondroitin polymerizing factor 2 (CHPF2), 


XM_003214253.1 


Meleagris gallopavo 


2.00E-178 


Reticulon 4 (RTN4) 


XM_003640893.1 


Gallus gallus 


3.00E-143 



Accession numbers are from the National Center for Biotechnology Information (NCBI) nonredundant nucleotide database. 
^Expected (E) value is a parameter when searching a sequence of particular size in a database (http://www.ncbi.nlm.nih.gov/BLAST/tutorial/Altschul- 
Lhtml). It is used to test the hypothesis of a random match to the bases in the database, the closer the E-value to zero, the lower the possibihty of 
mistake in a gene identity. 



and Ectonucleotide pyrophosphatase/phosphodiesterase 2 
(ENPP2). 

Functional annotation 

Based on the gene ontology (GO) tems, these genes 
were functionally annotated. As illustrated in Figure 1, 
majority of the biological processes (such as anatomical 
structure development, biosynthetic process, cellular 
nitrogen compound metabolic process, transport, cell 
differentiation, cellular protein modification process, signal 
transduction, small molecule metabolic process, anatomical 
structure formation involved in morphogenesis, 
transmembrane transport, cell death, translation, and DNA 
metabolic process), some cellular components (such as cell, 
intracellular, organelle, cytoplasm, mitochondrion, nucleus, 
protein complex and plasma membrane), and some 
molecular functions (such as kinase activity, ion binding 
and enzyme regulator activity) were presented in both 
libraries. Other GO terms, such as transmembrane transport. 



ribosome biogenesis, growth, endosome, nucleolus, ligase 
activity, lyase activity, lipid binding, peptidase activity and 
methyltransferase activity were identified only in the 
forward library. Vice versa, some GO terms were identified 
only in the reverse Ubrary, such as cytoskeleton 
organization, cell motility, lipid metabolic process, cell 
morphogenesis, cell cycle, protein targeting, membrane 
organization, neurological system process, locomotion, 
developmental maturation, vesicle-mediated transport, 
embryo development, cell prohferation, cell-cell signaling, 
plasma membrane organization, immune system process, 
proteinaceous extracellular matrix, cytosol, Golgi apparatus, 
endoplasmic reticulum, extracellular region, transferase 
activity, transferring acyl groups, nucleic acid binding 
transcription factor activity and cytoskeletal protein binding. 

Confirmation of gene expression witli quantitative 
RT-PCR 

To verify the results of SSH, we analyzed the expression 
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Table 3. Identities of differentially expressed cDNA clones from the SSH reverse library 



NCBI BLAST gene name 


Accession number' 


Species 


E-value^ 


IVTiisrlpVilinH-liVp f^T^ro^ionViilfi^ fA/IRNTT \\ frutrirrint vnrinnt 1 


NM_ 


_001 163338.1 


IrHllll'k (Tnllll'k 
vJ till Li L> ^^llLtL> 


0 ooF-i-nn 




XM_00 1234600.2 


UailUs gall Us 


1 OOF 1 A9. 

1 .UUl-i IfO 


Nralnha^-arfttvlfransfprasp 90 N^atT? ratalvtir siiliiinit ^^N^AA9fl^ 


XM_422177.3 


Oalliis (Tfllliis 

VJUJ-lUa gOJ-lUa 


5 OOE-57 


^nlntp parripr "ffimilv 1 ( crliiil hioli affinitv ffliitamatp trnn^nnrtpr^ mpmhpr (^J f^l A'^'t 

OL'lllLt' y^tU.Ll\/L Idlllll^ 1 yglldl lllgll alllllliy glUUCUlldLt' lLalliilJ\JH\>L lll\/llUJ\>L J Li^J J 


XM_425011.2 


vJdllLJS gullLJS 


0 


l^L.'l^JllU^lCVJtlUC lJYl\J|JllLJSlJlla.LtlSC/ UilUaUllUUJ-C&lClCl&C ^ \LIii.yi:r^.f 


NM_001 198662.1 




Q noF-i 10 

o .vjyc 1 1 o 


NnplpfiT" fVirtor T/R ('NPTR^ fraiTirriTit vnrifint 9 


NM^ 


_001 190738.1 


l-Tntnn ^aniPiiQ 


0 


vniirhpr TPIVTR 71^7 NAOl-r HphvHrnfrprifl^ip <:nhiinit 9 rnH9^ stptip initnrhnnHrial 


EU585683.1 


Anser rossii 


2.00E-115 




NM_ 


.205086.1 


VJaiiUa gdiiua 


l.OOE-34 




NM_ 


.001030971.1 


Gtillus gtillus 




iVllClUiiUllllaI-aasUl.-iaLCU piUlCill J \iyirt\ij) 


NM. 


.001012784.2 


ItOI ill C fTollll c 

kJailUs gclllUa 


A 


TntraflacTpllnr tmnQnoTt 74. linmnlrxT ('t~'lilamvHnmf\naQ'^ /" i h' 1 "74 ^ 


XM. 


.003643070.1 


Oalliid (Tfilliid 


1 OOF-50 


Trsnsmcinbraric cnip24 domain trsfficlciiig protein 2 (TMED2) 


NM_ 


.001006186.1 


GbIIus gsllus 


2 OOE-76 


Prnfilin 2 CPFN91 


NM_ 


.001079760.1 


VJ ail LIS gullLtS 


0 


Amyloid bctti precursor- like protein 2 (APLP2) 


NM_ 


J)01006317.2 


G^llus gtillus 




V^ydlll UCpCIlUCllL JsJ-llabC J, IC^UiaLUiy aUUUIlIl 1 IJJJ J ) L-* Ix 1 ^ 


NM. 


.003885.2 


nuillU sapiCIlb 


l.OOE-49 


Rihnsomal nrotein T '^9 rRPT ^Q"* 


NM. 


.204272.1 


rrallllS (TflllllS 


0 00000005 


Vacuolar protein sorting 37 iiomolog B 


XM_ 


.415126 


Gallus gallus 


4.00E-103 


Vesicle-associated membrane protein 4 (VAMP4) 


XM. 


.001233851.2 


Gallus gallus 


7.00E-21 


Rho guanine nucleotide exchange factor (GEF) 37 (ARHGEF37) 


XM. 


.414480.3 


Gallus gallus 


l.OOE-78 


Acetyl-CoA acetyltransferase 1 (ACATl) 


XM. 


.417162.3 


Gallus gallus 


4.00E-57 


Calmodulin binding transcription activator 1 (CAMTAl) 


XM_ 


.417530.3 


Gallus gallus 


0 


SPARC related modular calcium binding 2 (SMOC2) 


XM_ 


.419600.3 


Gallus gallus 


2.00E-123 


Mitochondrion, complete genome 


AF363031.1 


Anser albifrons 


6.00E-170 


hypothetical protein 


AJ720212.1 


Gallus gallus 


l.OOE-172 


Dynein, axonemal, light chain 1 (DNALl ) 


NM_ 


_001 199681.1 


Gallus gallus 


3.00E-158 



' Accession numbers are from the National Center for Biotechnology Information (NCBI) nonredundant nucleotide database. 

^ Expected (E) value is a parameter when searching a sequence of particular size in a database (http://www.ncbi.nlm.nih.gov/BLAST/utorial/Altschul- 
Lhtml). It is used to test the hypothesis of a random match to the bases in the database, the closer the E-value to zero, the lower the possibiUty of 
mistake in a gene identity. 



of 11 ESTs by qRT-PCR. The fold difference between the 
laying group and the ceased group was calculated by 2"'^'"' 
method. As illustrated in Figure 2, the qRT-PCR results for 
the most part, consistent with the SSH results. For the 6 up- 
regulated genes, 4 of them were over-expressed in the 
laying period, whereas 4 out of 5 down-regulated genes 
were under-expressed. Particularly, SYTl, STMN2, ENPP2, 
and CDK5R1 exhibited significantly differentially 
expressed in laying period pituitary compared with that in 
ceased period pituitary. 

DISCUSSION 

In the present work, we employed SSH method to 
identify differentially expressed genes in the pituitary gland 
of Huoyan geese between laying period and ceased period. 
Some genes involved in signal transduction (such as 
CNTNAP2, STMN2, PAK2, and SMOC2), calcium ion 
binding (such as CAMTAl, PCDH9, CDK5R1, and SYTl), 
MAPK signaUng pathways (SPREDl), I-kappaB 
kinase/NF-kappaB signaling pathways (VAPA), and 
transport (such as VAMP4, SLC5A11, SLC4A4, COXl, 
and COX3) were identified. Subsequently, eleven of those 
including PAK2, PCDH9, SPREDl, STMN2, SYTl, VAPA, 



CAMTAl, ENPP2, SMOC2, VAMP4, and CDK5R1 were 
selected for further validating using qRT-PCR. The results 
were mostly consistent with the SSH results. Interestingly, 
several of these genes might be involved secretion of 
hormones. 

One of the major findings was that several genes related 
to exocytosis were found differentially expressed in our 
study. Exocytosis is a key biological process that controls 
the neurotransmission and release of secretory products 
from neurones and other secretory cell types. 
Neurotransmitters, hormones, or other secretory products 
are packed in vesicles, a number of these vesicles fuse with 
the surface membrane, during both nonactivated and 
activated phases, to release secretory product into the 
extracellular space (Knight and Baker, 1987). 
Neurotransmitters and hormones release is triggered by 
Ca^"^ binding to a presynaptic Ca^"^ sensor that induces 
synaptic vesicle exocytosis with a high degree of Ca^"^ 
cooperativity. Synaptotagmin-1 (SYTl) has been identified 
as the primary Ca^^ sensors for synaptic vesicle exocytosis 
and a major transducer of Ca^"^ signaling in membrane 
fusion events and regulated secretion (Koh and Bellen, 
2003; Kreft et al., 2003; Shin et al., 2009). There is 
evidence that SYTl plays a physiological role in secretion 
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Figure 1. Gene ontology (GO) classification of the differentially expressed ESTs in the forward and reverse SSH Ubraries. Based on the 
GO terms, the ESTs are categorized into putative functional groups. The black and white columns represent the up- and down-regulated 
genes from forward and reverse libraries, respectively. 



by differentiated pituitary cells. Genetic experiments in 
mice have demonstrated that SYTl mutants have defects in 
regulated secretion (Koh and Bellen, 2003). Mice with 
homozygous disruption of the SYTl die shortly after birth 
and have defects in neurotransmitter release from 
hippocampal neurons (Geppert et al., 1994). The pituitary- 
specific transcription factor (POUIFI) is a factor that bound 
to and activated the PRL and GH promoters, and is also 
important for proper development of the pituitary cells that 
express PRL, GH, and TSH. Study provides evidence that it 
can bind to a specific site in the SYTl gene and that this 
binding contributes to activation of SYTl expression. It 
seems likely that activation of SYTl gene expression is part 
of a mechanism mediating POUlF-induced differentiation 
of pituitary cells and presumably contributes to the 
endocrine/secretory phenotype (Howard et al., 2009). 
Production and secretion of hormones by the pituitary 
involve highly orchestrated intracellular transport and 
sorting steps. Peptide hormones destined for regulated 
secretion are packaged into secretory granules that bud from 
the frans-Golgi network (TGN). VAMP4 (Vesicle- 
associated membrane protein 4), a member of VAMP 
proteins implicated in vesicle docking with membranes, has 



been suggested to play a role in mediating trans-Go\g\ 
network (TGN) trafficking. VAMP4 marks a novel 
vesicular sorting pathway during granule membrane 
remodeling, and Synaptotagmin TV follows the trafficking 
patterns of VAMP4 and modulates Ca^"^-triggered 
exocytosis during granule maturation (Steegmaier et al., 
1999). 

It is well known that Ca^"^ is the primary intracellular 
signaling molecule controlling the fusion of secretory 
vesicles with the plasma membrane to release transinitters 
from neurons and hormones from endocrine cells (Martin, 
2003). Increase in intracellular Ca^"^ concentration is a key 
signal to trigger exocytosis in most neuronal and endocrine 
tissues. Pituitary hormone secretion is a Ca^*-dependent 
process. Depriving extracellular Ca^"^ by incubation in 
media containing Ca^"^ buffers such as EGTA inhibits 
GnRH-induced LH secretion (BATES and CONN, 1984), 
GRF-induced GH secretion (Sheppard et al., 1980) and 
TRH-induced PRL secretion (Thorner et al., 1980). 
Furthermore, this inhibition can be reversed by elevation of 
extracellular calcium. In this study, except SYTl, severed 
other genes which involved calcium ion binding were 
identified differentially expressed. Protocadherins (PCDHs) 
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Figure 2. Quantitative real-time PCR validation. The gene expression level of 11 randomly selected ESTs, including 6 from SSH forward 
library (the black column) and 5 from SSH reverse library (the white coluiim), was tested using qRT-PCR. The relative gene expression 
level in the laying period was normalized to the ceased period. Data are presented as mean+SE. Asterisks denote significant gene 
expression differences between the laying period and ceased period, as determined by the student's t tests (* p<0.05). 



constitute the largest subgroup within the cadherin family of 
calcium-dependent cell-cell adhesion molecules, and have 
been suggested to play role in the formation and 
maintenance of the synaptic connections (Takeichi and Abe, 
2005). Study demonstrated some PCDHs were involved in 
regulation of LHRH neuronal connections (Varju et al., 
2009). Calmodulin can trigger exocytosis of already primed 
vesicles in endocrine cells; it also binds to the Ca^"^ sensor 
synaptotagmin, and could be involved in synaptotagmin- 
related events (Langley and Grant, 1997). CalmoduUn 
binding transcription activator 1 (CAMTAl) was 
differentially expressed in our study, it may be involved 
alterations of intracellular Ca^"^ concentration and play a 
major role in anterior pituitary cell stimulus-secretion 
coupling. 

Remarkably, the expression of Stathmin-2 (STMN2) 
was identified significantly up-regulated in laying period. 
Stathmin is the founding member of a family of proteins 
that play critically important roles in the regulation of the 
microtubule cytoskeleton, its expression and 
phosphorylation are regulated throughout development and 
in response to extracellular signals regulating cell 
proUferation, differentiation and functions (Rubin and 
Atweh, 2004). It is highly expressed in the nervous and 



neuroendocrine systems and its phosphorylation could be 
related to regulation of cell growth arrest, terminal 
differentiation and hormone secretion (Rowlands et al., 
1995; Drouva et al., 1998). Studies demonstrated that 
Stathmin play a potential role in the regulation of hormone 
secretion in rodent pituitary and insulinoma cell lines 
(Rubin and Atweh, 2004). Interestingly, the mRNA 
expression level of STMN2 was found to be significantly 
higher in the hypothalamus and pituitary gland of high egg 
production chicken strains (Shiue et al., 2006). Its role in 
regulation of laying egg needs to be investigated in the 
future. 

Furthermore, several genes involved in regulation of 
kinase activity were identified differentially expressed in 
this study. Sprouty-related proteins (SPREDs) with EVHl 
(enabled/vasodilator-stimulated phosphoprotein homology 
1) domain are membrane-associated suppressors of MAPK 
signaling, a major regulator of cellular proliferation and 
differentiation. SPRED2 was identified as a negative 
regulator of the Hypothalamic -Pituitary- Adrenal Axis, 
SPRED2 deficiency leads to an increased MAPK signaling, 
which results in an augmented its secretagogue 
corticotrophin releasing hormone (CRH) promoter activity, 
the subsequent CRH overproduction causes an up- 
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regulation of downstream hypothalamic -pituitary-adrenal 
(HPA) hormone secretion (Ulkich et al., 2011). SPREDl 
has been shown to be expressed predominantly in the brain 
and enriched in the central nervous system, where it 
mediates cortical development, neural stem cell 
proliferation, and vesicular trafficking (Phoenix and Temple, 
2010). In this study, the differential expression of SPREDl 
may be associated with its vesicular trafficking and negative 
regulation of the Ras-MAPK-ERK pathway during 
production and secretion of hormones in the pituitary gland. 
CDK5R1, also named p35 was the first neuronal-specific 
regulatory subunit of CDK5 (Cyclin-dependent kinase) 
identified via its physical interaction with CDK5, activating 
CDK5 after direct binding (Ko et al., 2001). Study suggest 
that CDK5/p35 play critical roles in neural migration, 
differentiation, and synaptogenesis in immature neurons. 
Cdk5/p35 also was indicated inhibits neurotransmitter 
release through the phosphorylation of P/Q-type voltage- 
dependent Ca^'^ channel and down-regulation of the channel 
activity (Tomizawa et al., 2002). The brain 
neurotransmitters involve the mechanism triggering the 
release of the hypothalamic hormones. Changes in 
neurotransmitter output and in particular alterations in 
secretion of the monoamines, dopamine, noradrenaline and 
serotonin, have been associated with many hormonal 
changes in mammals (Chen et al., 1974; Advis et al., 1978). 
The effects of neurotransmitters on release of pituitary 
hormones in birds have been studied. Monoamines were 
confirmed to be capable of affecting pituitary release of 
PRL and GH in a pituitary-hypothalamus in pigeon (Hall, 
1982). 

In summary, we have shown gene expression profiles in 
the pituitary gland of laying period and ceased period 
Huoyan geese. Particularly, SYTl and STMN2 which play 
an important role in release of pituitary hormones were first 
to demonstrate substantially over-expressed in the laying 
goose pituitary gland. Further demonstration of the 
functions of genes will add great value to the imderstanding 
of the reproductive biology of the goose. 
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